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[57] ABSTRACT 

A polarization filter can maximize the signal-to-noise 
ratio of a polarimetric SAR and help discriminate be- 
tween targets or enhance image features, e.g., enhance 
contract between different types of target. The method 
disclosed is based on the Stokes matrix/Stokes vector 
representation, so the targets of interest can be extended 
targets, and the method can also be applied to the case 
of bistatic polarimetric radars. 

7 Claims, 5 Drawing Sheets 
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METHOD FOR PROVIDING A POLARIZATION 
FILTER FOR PROCESSING SYNTHETIC 
APERTURE RADAR IMAGE DATA 

5 

ORIGIN OF THE INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected not to retain 10 
title. 

TECHNICAL FIELD 

The invention relates to a method for providing a 
polarization filter for processing synthetic aperture 15 
radar (SAR) image data, and more particularly for a 
filter for maximization of signal-to-noise ratio and a 
filter for discrimination between areas of different tar- 
get types by utilizing the vector nature of the return 
waves in order that information about the surface con- 20 
tained in the polarization properties of the scattered 
return waves may be recovered. 

BACKGROUND ART 

Conventional imaging radars operate with a single 25 
fixed-polarization antenna for the radio frequency sig- 
nal transmission and reception. As a consequence, for 
every resolution element (pixel) in the image, the scat- 
tered wave (a vector quantity) is measured as a scalar 
quantity and any additional information about the sur- 30 
face contained in the polarization properties of the re- 
flected signal is lost. To insure that all the information in 
the reflected signal is retained, the reflected signal must 
be measured as a vector, i.e., both the amplitude and the 
phase should be measured. The greater information 35 
derived from the polarized measurements provides a 
more complete description of the backscatter phenom- 
ena of the target area. This greater information can, for 
example, be used in target discrimination, target classifi- 
cation or feature enhancement. This invention relates to 40 
techniques for utilizing polarization information to en- 
hance certain characteristics in SAR images. 

Determination of the optimal polarization state to 
enhance an image has been extensively studied in the 
past. The scattering matrix co-polarization and cross- 45 
polarization nulls represent a solution when the returns 
from a single point target need to be minimized. A. B. 
Kostinski and W. M. Boemer, “On the Polarimetric 
Contrast Optimization,” IEEE Trans. Antennas and 
Propagation, Vol. AP-35, No. 8, pp. 988-991, August 50 
1987, based their analysis on the Graves power matrix 
to determine the optimum transmit and receive polar- 
izations. This technique is only applicable to maximize 
the contrast between two specified point targets and has 
been applied to polarimetric radar images by estimating 55 
an equivalent scattering matrix representation for an 
extended area in an image. B. James, A. B. Kostinski 
and W. M. Boemer, “Polarimetric Matched Filter for 
POLSAR Image Interpretation of Ocean Surface Scat- 
ter,” Proc. IGARSS ’88, P. 67, Edinburgh, United 60 
Kingdom. 

The accuracy of this technique for extended targets is 
unknown since representation requires either an aver- 
age Stokes matrix [J. J. van Zyl, A. Zebker and C. 
Elachi, “Imaging radar polarization signatures: Theory 65 
and observation,” Radio Science, 22(4), pp. 529-543, 
July/August 1987; and J. J. van Zyl, C. H. Papas and C. 
Elachi, “On the optimum polarizations of incoherently 


2 

reflected waves,” IEEE Trans, on Antennas and Propa- 
gation,” Vol. AP-35, No. 7, July 1987] or an average 
covariance matrix [J. A. Kong, A. A. Swartz, H. A. 
Yueh, L. M. Novak and R. T. Shin, “Identification of 
Terrain Cover Using the Optimum Polarimeter Classi- 
fier,” J. Electromagnetic Waves and Applications, Vol. 
2, No. 2, pp. 171-194, 1988]. Both these representations 
use the second order statistics of the scattering matrix. 
G. A. Ioannidis and D. E. Hammers, “Optimum An- 
tenna Polarizations for Target Discrimination in Clut- 
ter,” IEEE Trans. Antennas Propagat., Vol. AP-27, 
No. 3, May 1979, introduced a method based on La- 
grangian multipliers to solve for the optimal polariza- 
tion using the Stokes matrix. However, for some cases 
their solutions violate the constraint that the Stokes 
vector for the receive antenna must be fully polarized. 
Additionally, A. A. Swartz, H. A. Yueh, J. A. Kong, L. 
M. Novak and R. T. Shin, “Optimal Polarizations for 
achieving Maximum Contrast in Radar Images,” J. 
Geophys. Res., Vol. 93, No. B12, pp. 15252-15260, 
December 1988, developed a parallel method based on 
the covariance matrix. However their analysis is re- 
stricted to the backscatter case. 

STATEMENT OF THE INVENTION 

An object of this invention is to provide a polariza- 
tion filter for maximization of the signal-to-noise ratio of 
a single target return or maximization contrast between 
two target types, particularly between man made urban 
areas and natural nonurban areas (forest, grass, and 
ocean). 

In accordance with the present invention, polariza- 
tion filters based on Stokes matrix/Stokes vector repre- 
sentation maximizing the signal-to-noise ratio are devel- 
oped for different noise characteristics and maximizing 
the contrast between different target types by maximi- 
zation of the power ratio between target types. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. la illustrates the total SAR power image (also 
known as the Span of the Stokes matrix). 

FIG. lb illustrates the outlines of different areas of 
target types and the optimum filter image of the SAR 
power image of FIG. la. 

FIG. lc illustrates the enhancement factor image 
which is the ratio of the two other images. 

FIG. Id shows the outlines of different areas of target 
types and gives the position of two studied sites; No. 1 
is a forested area site and No. 2 is an urban area site. 

FIG. 2 illustrates an optimization signature maximiza- 
tion of signal-to-noise ratio for an urban area. 

FIG. 3 illustrates an optimization signature maximiza- 
tion of signal-to-noise ratio for a forested area. 

FIG. 4 illustrates an optimization signature maximiza- 
tion of contrast ratio between urban and forested areas. 

FIG. 5 is a graph of enhancement factor variation 
with range line. The line 0 is the near range and line 
1000 is the far range of the SAR image. 

FIG. 6 is a graph of enhancement factor variation 
with size of training area. The training area contains n 2 
pixels. 

FIG. la illustrates the total power image (also known 
as Span of the Stokes matrix) of the Blackwater River 
area. 

FIG. lb is the optimum filter image of FIG. la , and 

FIG. lc is the enhancement factor image of FIG. la, 
which is the ratio of the two other images. 
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FIG. 8 illustrates an optimization signature maximiza- 
tion of signal-to-noise ratio for an urban area with sym- 
metrized noise. 

DETAILED DESCRIPTION OF THE 5 

INVENTION 

A systematic theoretical analysis and development of 
a polarization filter is described for two applications: (1) 
maximization of signal-to-noise ratio and (2) discrimina- 
tion between two target types. The method employed 10 
finds the optimum receive antenna polarization analyti- 
cally, but relies on a numerical solution to find the opti- 
mum transmit antenna polarization. The analysis uses 
the Stokes matrix representation, and therefore can be 15 
employed to analyze the partially polarized scattered 
field from extended targets and can also be applied to 
the bistatic case. 

Examples of the technique are presented for the 
monostatic NASA CV990 polarimetric L-band radar 2 o 
which alternately transmits illuminating waves with 
horizontal (h) and vertical (v) polarization while the 
receiver employs two channels, one for the horizontally 
polarized scattered waves (h') and the other for the 
vertically polarized scattered waves (v') in order that 25 
the data received and stored may be processed in a 
bistatic scattering matrix 


N = 


A k’h N hv 
Ny'h N v ' v 


( 1 ) 


Each element of N is a complex random variable with 
zero mean and a known variance cr 2 (the noise power). 
Assume that the elements of M are uncorrelated and 
that the noise of each radar channel has the same power. 
The mathematical formulation is as follows: 


<NijNkl> = cr 2 if i = k and j — 1 
= 0 otherwise 

<j V, y> = 0 for all i and j 


where < > denotes the spatial average. 

The corresponding average Stokes matrix is of the 
form: 


M N = 


cr 2 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 


(2) 


Sh'h 5/,-v 

Sy’h Sv'v 


30 


where: Skh is the complex ratio of the electric field of 
the horizontally polarized part of the scattered wave 35 
(h') and the horizontally polarized part of the illuminat- 
ing wave (h); 

Sh'v is the complex ratio of the electric field of the 
horizontally polarized part of the scattered wave (h') 
and the vertically polarized part of the illuminating 40 
wave (v); 

S v 'h is the complex ratio of the electric field of the 
vertically polarized part of the scattered wave (v') and 
the horizontally polarized part of the illuminating wave 
(h); and 45 

S v 'y is the complex ratio of the electric field of the 
vertically polarized part of the scattered wave (v') and 
the vertically polarized part of the illuminating wave 

(v). 

The terminology and definitions used here and here- 
inafter have previously been defined by J. J. van Zyl, H. 

A. Zebker and C. Elachi, (July/August 1987) and are 
incorporated herein by reference. 

Image enhancement filters maximizing the signal-to- 
noise ratio are developed for different noise characteris- 
tics and different target types (urban and forest). A filter 
is also developed to maximize the power ratio between 
urban and natural (forest or ocean) targets for optimum 
discrimination between the types of targets. Results are 
presented to demonstrate that the filter maximizing the 
contrast between urban and ocean areas is essentially 
the same as the one maximizing the contrast between 
urban and forested areas. 


Matched Filter in Presence of Background Noise 

Let N be a scattering matrix characteristic of the 
received noise. 


Equation (2) states that the average noise signal is unpo- 
larized when the noise is uncorrelated. Therefore, the 
noise merely contributes a constant power level to the 
signal [van Zyl, Papas and Elachi (July 1987)]. 

To find the polarization configuration of the receive 
and transmit antenna maximizing the signal-to-noise 
ratio (SNR), let A be the average Stokes matrix charac- 
terizing a target of interest, a M^ybe the average Stokes 
matrix for the noise. Also let S r and S t be the Stokes 
vectors characterizing some arbitrary polarization con- 
figuration of the receive and transmit antennas. 


Siq * S/o 
5/ — and S r = 

s { s r 


where (3) 



Sri 


St\ 

S r = 

Sr2 

and St — 

Sa 

50 

Sr 3 


St 3 


(3) 


Because both Stokes vectors represent an antenna con- 
figuration, a necessary condition [van Zyl, Papas and 
Elachi, (July 1987)] is: 


55 


Sfi = • 


/ and si) = . 


(4) 


It can be shown [van Zyl, Papas and Elachi, (July 
1987)] that the signal-to-noise ratio SNR is given by: 


SNR = t 


S r T A St 


(5) 


V M N S t 


65 Without loss of generality, it can be assumed that 
S r g=Srf)= 1 since these terms cancel in the ratio. 

The received power for a given antenna polarization 
can be maximized or minimized by maximizing or mini- 



5,057,843 


5 

mizing S and S t0 • However, in practice S/o and S/o are 
proportional to the gain of the antennas, which for a 
given wavelength is determined by their sizes and can- 
not easily be varied. Consequently, the present inven- 
tion uses signal polarization characteristics to extract 5 
information from SAR data assuming that the gain of 
the transmit and receive antennas are fixed. Therefore, 
Equation (4) can be rewritten as: 


s r s r = 1 and s t s t = I (6) 1Q 

Assuming uncorrelated noise, the denominator of 
Equation (5) is constant and equal to cr 2 . Therefore, 
maximizing the signal-to-noise ratio is equivalent to 
maximizing the numerator. To optimize SNR, the aver- ^ 
age Stokes matrix A is written as follows: 



where u denotes the transpose of u and u, v are 3-ele- 
ment real vectors and Q is a 3x3 real matrix. 

Using Equation (7) for A, an alternate form for the 
numerator of Equation (5) is: 

S r T A S t — m -f u • s t + v • s r + s r * Qs t 



This is the expression to be maximized under the con- 
straints expressed in Equation (6), namely the assump- 
tion that the gain of the transmit and receive antennas 
are fixed. However, the solution of constrained maxima 
for any one of the variables may be extremely difficult. 

It is therefore desirable to use the Lagrangian multiplier 
method to solve this problem. After the necessary dif- 
ferentiation, the s/and s r that maximizes the SNR is the 35 
solution of the equations: 


u + Qs r = M s t (9) 

v + Qs t = (10) 

where Ai and X2 are the Lagrangian multipliers. 

If the transmit antenna is known (i.e., s r is given), the 
receive antenna polarization which maximizes the SNR 
is determined from Equations (6) and (9) to be as fol- 
lows: 


v+ Qs t (ll) 

p+&,h 

where || || denotes the norm. Physically, Equation (11) 50 
states that the power is maximized when the polariza- 
tion of the receive antenna is matched to the polariza- 
tion of the scattered wave incident upon the receive 
antenna. 55 

The transmit antenna polarization Sfinust then satisfy: 

(Q Q-M*ifrt=-Q v-X 2 «. (12) 

When Equation (1 1) is satisfied, the numerator of Equa- 
tion (5) can then be written: 60 

P max— mJ r u-s t +\UQst\ (13) 

where F ma x is the maximum obtainable scattered power 
given the transmit antenna polarization. 65 

It is hard to analytically solve for the transmit polar- 
ization which would maximize V max . Instead, this opti- 
mum transmit polarization may be found numerically. 


6 

as will be shown below under the heading Numerical 
Examples. 

Matched Filter as a Discriminator Between Two 
Target Classes 

One objective of this invention is to find the polariza- 
tion of the receive and transmit antennas which maxi- 
mizes the ratio of signal power scattered by one type of 
target to that scattered by another type of target. Let 
F(D represent the Stokes matrix characteristic of target 
No. 1 and F( 2 ) represent the Stokes matrix characteristic 
of target^No. 2^The problem then is to find the Stokes 
vectors S r and S t representing the receive and transmit 
antenna polarization realizations which maximize: 

£ i) _ s r T l*»St (14) 

/* 2 ) “ S r T i* 2 > S t 

But first it is necessary to solve the following problem. 
Assuming the transmit antenna polarization is fixed, 
only the receive antenna is to be operated upon to maxi- 
mize the power ratio given by Equation (14). Si and S2 
are defined as follows: 


S { = jtfDS/ = 


S2 = JF^^St — 


I'Soi I 

I s \ I 

|S°2 I 

I s 2 I 


(15) 


(16) 


S r , the Stokes vector of the receive antenna, can be 
written as: 



(17) 


As before, since S r represents the Stokes matrix of an 
antenna, it is a completely polarized wave and hence 
must satisfy: 

s 2 r0 = i r . ir =l (18) 

It is clear that to maximize the ratio given by Equation 
(14), it is necessary to maximize the contrast 


„ x *Soi + *i (I 9 ) 

C(Sr) ~S02 + Sr- *2 

Since C(s r ) is the ratio of two received powers, it is 
reasonable to insist that: 


C(s r )=0. (20) 

To maximize the contrast between the two targets, one 
has to maximize the return from the one target, while at 
the same time minimizing the return from the other 
target. For the moment, assume that there is no antenna 
polarization which could cause the denominator of 
Equation (19) to become zero and then generalize the 
result. To optimize C(s r ) subject to the constraint condi- 
tion given by Equation (18), Lagrangian multipliers are 
introduced. After performing the necessary differentia- 
tions, it is found that the optimum polarizations are the 
solutions to 
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-continued 


/ \ ' ' 
( Sq] + s r * s\ ) 

51 - I2 Uo 2 + ^-2 

where /x is the Lagrangian multiplier. Introducing a 
shorthand term 

_ 5qi + s r • s\ (22) 

“ 502 + Sr - s 2 ’ 


Equation (21) can be rewritten to read 

s(a)—iis r 
s(a)^si — aS2- 


( S 01 S 02 — S\ ■ s 2 

2 

s oi - s \ ■ s \ 

^ S 52 - S2-S2 

«Sq2 — 52 - 52 


From Equations (19) and (20), it is clear that a is also 
the ratio to optimize, and therefore Equations (28 a and 
JO b) also define the optimum values of the contrast ratio, 
C(s r ). It can be proven that the contrast ratios given in 
Equations (28) are real, positive and definite values. 

Substituting Equations (28a) and (286) in Equation 
(23) and substituting the resulting expression in Equa- 
ls tion (25), it is found that only the solutions with positive 
signs lead to self-consistent results for the contrast ratio. 
Therefore, the optimum antenna polarizations are: 


From Equations (23) and (18) it follows that 


±v 


5(a) • s(a ) 


20 

(24) 


■Sterner*) 

rmax | | .y(a/7jcjt) I j 

and 


(29a) 


and the optimum polarizations are therefore described 
by 

*-±., M . < 25 > 

N *a) • j(a) 

Notice that s r , as expressed in Equation (25), always 
represents a fully polarized wave consistent with the 
constraint in Equation (18). A previous analysis by Ion- 
nadis and Hammers (1979) found the optimum receive 
antenna Stokes vector: 


sia min ) (29b) 

hmin - | |i(a m/B )| | ' 

Now notice that when So2 2 =r=S2-S2, the scattered wave 
in the denominator of Equation (19) is completely po- 
larized. The maximum contrast in this case is obviously 
obtained when an antenna polarization which nulls the 
received power from the completely polarized wave is 
used to receive the scattered waves, leading to an infi- 
nite contrast ratio. The receive antenna polarization is 
thus chosen to be orthogonal to the completely polar- 
ized scattered wave: 


25 


30 


35 


Sr = 4- JKJt'f - ^ 2) )S| = 4- S 1 As, 


(26) 


where p is a Lagrangian multiplier. In general, Si and ^ 
S2 represent partially polarized waves if FO) and F( 2 ) are 
Stokes matrices of extended targets. There is therefore 
no guarantee that the difference S r would be a fully 
polarized Stokes vector. 

In Equation (25), a is still unknown. That unknown 45 
can be solved by substituting the expression for s r in 
Equation (25) into the definition of a given in Equation 
(22). This produces the optimum values of a that are the 
roots of the quadratic: 

50 

(S 02 “ si * S2)a 2 - 2(5oiSo2 - *1 * s 2 )a + 0Sqi - s\ • s\) = 0 (27) 


If So2 2 9^S2 S2, the optimum values of a are 

55 

S01S02 ~ s\ -S2 ( 28 a) 

a max ~ ry + 

Sq 2 ~ S2-S2 


^ SoiSo2 “ S 1 * *2 ^ 


•Sn? — 


*2 • S2 


£01 - *1 ♦ *1 
Sq 2 - s 2- S 2 


60 


and 


a min 


SqiSqi ~ S\ • S 2 
*5q2 — S2 ' S2 


65 

(28b) 


s rmax 


s 2 

N 52 * S 2 


(30) 


To find the minimum contrast ratio for this case, return 
to Equation (23) and find that the quadratic reduces to 
the first order polynomial with 

Sq\-SI-51 (31) 

am ' n= 2 (SoiSo2 - *1 ’*2) 

Notice that a m /„ as given by Equation (31) is the ratio of 
two positive quantities and therefore satisfies Equation 
( 18 ). 

Note that if the scattered wave in the numerator of 
Equation (18) is completely polarized, the minimum 
contrast ratio is zero which is achieved by using the 
antenna polarization which nulls the received power 
from this completely polarized wave. When both waves 
are completely polarized, the minimum contrast (which 
is zero) is achieved by using the antenna polarization 
which nulls the received power from the completely 
polarized wave in the numerator of the contrast ratio. 
The maximum contrast (which is finite) is achieved by 
using the antenna polarization which nulls the received 
power from the completely polarized wave in the de- 
nominator of the contrast ratio. 

Finally, it should be noted that if the F( 2 ) matrix ap- 
pearing in the denominator of Equation (14) is a noise 
matrix as defined in Equation (2), then from Equations 
(16) and (21) s 2 =0 and 
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Si ( 32 ) 

i,_ Hull 

where s\ = v+Qs r . This result agrees with Equation (2). 

The equations derived in this section describe the 
optimum receive antenna polarization for maximum 
contrast ratio given the transmit antenna polarization. 
Furthermore, Equation (28a) expresses the maximum 
contrast ratio obtainable given the transmit antenna 
polarization state. The next logical step is to vary the 
transmit antenna polarization to maximize the contrast 
over all possible configurations. Attempts to find an 
analytical solution for the optimum transmit polariza- 
tion lead to extremely complicated nonlinear equations. 
Instead of solving these, this invention presents a sim- 
pler hybrid method in which the optimum transmit 
polarizations are found numerically by calculating the 
maximum and/or minimum contrast ratio. The corre- 
sponding receive polarizations are calculated using 
Equations (28) and (29) for all possible transmit polar- 
izations. This hybrid method is illustrated with several 
examples in the following section. 

Numerical Examples 

The procedures developed in the previous sections 
were tested with data acquired with the NASA CV990 
SAR operating at L band. This radar is monostatic and 
emits alternately horizontally and vertically polarized 
waves and receives both returns with two co-located 
antennas, one horizontally polarized and the other ver- 
tically polarized. The system is therefore able to mea- 
sure a full scattering matrix for each of the 12x4 meter 
resolution cells. 

The first test was to minimize the impact of back- 
ground noise using the matched filter approach in the 
presence of background noise described above. A sec- 
ond test was to demonstrate the use of a polarimetric 
matched filter to discriminate between man-made tar- 
gets and natural targets. A SAR image shown in FIG. 1 
of San Francisco acquired in 1985 was used in both 
examples. It was chosen because of the diversity of the 
target types as well as the availability of the ground 
truth. FIG. la is the total power image, i.e., the (1,1) 
element of the Stokes matrix. FIG. lb is the matched 
filter output simulating 45° linearly cross polarized ra- 
dar. FIG. lc is the enhancement factor image given by 
Equation (34) below. FIG. Id is an outline of different 
major target areas, namely urban, forested and ocean. 
The data presented in the following examples have been 
relatively calibrated as that term is defined by Zebker, 
et al., (1986) “Imaging Radar Polarimetry from Wave 
Synthesis,*’ J. Geophys., Res., Vol. 92, No. Bl, pp. 
683-701, January 1987, but not absolutely calibrated. 

EXAMPLE A: SIGNAL-TO-NOISE RATIO 
MAXIMIZATION (VEGETATION, MAN-MADE 
TARGETS AND OPEN WATER) 

First the matched filter approach is applied to maxi- 
mize the signal to noise ratio. For a given target, the 
optimum receive polarization can be analytically com- 
puted for each transmit polarization. If implemented 
with sufficiently small increments, a numerical search 
over all polarizations of the transmit antenna will yield 
a global maximum of the signal-to-noise ratio as well as 
the corresponding receive antenna polarization configu- 
ration. This search process can be visualized through an 
“optimization signature” which is here introduced and 


10 

defined as a 3D plot where the two horizontal axes 
describe the polarization of the transmit wave, one axis 
being the ellipticity angle, the other the orientation 
angle. The maximum signal-to-noise ratio, calculated 
5 using Equation (12) at each transmit polarization is 
displayed along the vertical axis. The highest point of 
the surface will therefore be the global maximum of the 
signal-to-noise ratio while its horizontal coordinates 
describe the polarization of the optimal transmit an- 
tenna. The receive antenna polarization is computed 
with Equation (9). 

An optimization signature is displayed for two differ- 
ent types of targets. For the man-made target area No. 
15 2, shown in FIG. Id, the Stokes matrix is given by: 


291.94 17.35 29.42 14.66 

17.35 229.70 83.09 - 12.73 

29.42 83.09 - 42.83 - 22.89 

14.66 - 12.73 - 22.89 105.07 

The corresponding optimization signature indicates that 
the maximum signal-to-noise ratio occurs for a transmit 
25 polarization of about 10° orientation angle and 0° ellip- 
ticity angle. The exact numbers are \J/,= 10.0°; Xr=0.0°; 
tr=9.8 0 ;xr=-5 o . 

For the natural target area, the Stokes matrix is: 
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F</) = 


63.63 - 6.83 4.87 1.29 

- 6.83 32.17 1.64 - 3.82 

4.87 1.64 1.92 3.33 

1.29 - 3.82 3.33 12.29 


( 34 ) 


35 The corresponding optimization signature shown in 
FIG. 3 indicates that the maximum signal-to-noise ratio 
occurs for a transmit polarization of about 85° orienta- 
tion angle and 5° ellipticity angle. The exact numbers in 
^ this case are i(/,=85.0°; Xr=5.0°; i|/ r =84.7°; Xr= — 5.6°. 

EXAMPLE B: DISCRIMINATION BETWEEN 
MAN-MADE AND NATURAL AREAS. 

The matched filter approach described above as a 
45 discriminator between two classes is used as follows: 
Two typical target areas are chosen as training sites, the 
first a forested area No. 1 and the second an urban area 
No. 2 indicated in FIG. Id by small rectangles. The two 
areas are similar in both size and radar viewing angle. 
50 The corresponding polarization matched filter is com- 
puted using the two experimentally determined Stokes 
matrices. An optimization signature, as defined in the 
previous paragraph, is determined and presented in 
FIG. 4 . For this test, the vertical axis is the normalized 
55 contrast ratio defined by Equation (13) where the maxi- 
mum contrast ratio is normalized to 1. The peak of the 
surface occurs for a transmit polarization with 7.5° 
ellipticity angle and 47.5° orientation angle. This proce- 
dure is repeated for twenty target pairs. The polariza- 
tion of the optimal transmit and receive antennas is 
shown as a scatter plot in FIG. 5. The resulting orienta- 
tion and ellipticity angles cluster well around ^ = 45° 
and xr=0% i|/r=135 0 and xt= 0°. This corresponds to 
65 a linearly cross-polarized radar oriented at 45° from 
horizontal. As a measure of the performance of the 
filtering method, we define an enhancement factor as 
follows: 
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part S r T P* rb S t fj°{ (35) 

~Pf° r “ S r T Ff° r Si P u { \ b 

5 

The enhancement factor is therefore a measure of how 
radar polarization provides discrimination between tar- 
gets as a standalone tool (i.e., without utilizing the tar- 
get characteristics). 

FIG. 5 shows that for this example the enhancement io 
factor is about 3.0 dB and is range independent. The size 
of the training area does not seem to have an influence 
on the enhancement factor as shown in FIG. 6. This 
provides a qualitative representation of polarization 
image enhancement. 15 

This procedure was also applied to maximize the 
contrast between urban and ocean areas. The resulting 
optimum filter is found to be essentially the same as the 
one optimizing the contrast between urban and forested 
areas. FIG. 7 displays the effect of the filter on an ocean 2 q 
scene. FIG. la is the total power image, FIG. lb is the 
filtered image simulating a 45° linearly cross-polarized 
radar and FIG. 1c is the enhancement factor. This filter 
seems to be a good urban/natural discriminator inde- 
pendent of the type of natural targets. 25 

If the uncorrelated noise described with reference to 
Equations (1) through (12) is symmetrized before aver- 
aging, the corresponding Stokes matrix will have a 
different form. Let N be the symmetrized scattering 
matrix. We have: 3Q 

N hv =N vh =h(N hy +N yh ) (36) 


Therefore, the corresponding Stokes matrix, Mjv is 
written: 

35 


3a 2 0 


A 1 


0 

0 

0 


a 2 

0 

0 


0 

0 

a- 

0 


0 

0 

0 

a 2 


(37) 
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target is water, forest, or grass field. Results show that 
this optimum filter improves the contrast ratio by about 
3 dB, and that this enhancement is essentially indepen- 
dent of the incidence angle. 

The optimization signature discussed is a useful tool 
as it indicates both the variation and maxima of the 
contrast ratio as a function of the receive antenna polar- 
ization. The polarization resulting in the maximum con- 
trast ratio defines the optimum filter. The robustness of 
this filter is related to the variation of the optimization 
signature surface in the neighborhood of the maximum 
contrast polarization. For example, a rapidly varying 
surface implies that the filter performance degrades 
quickly when small errors occur in polarization realiza- 
tion. Such small errors can result from system imperfec- 
tions such as crosstalk. 

What is claimed is: 

1. In a polarimetric synthetic aperture radar imaging 
system, a method for polarization filtering in order to 
maximize the signal-to-noise ratio comprising the steps 
of 

assuming a constant gain of said transmit and receive 
antennas and assuming uncorrelated noise as a first 
step, 

then maximizing S r r A S?, as a second step where Sr 
and S r r are the Stokes vectors characterizing some 
arbitrary polarization configuration of said trans- 
mit and receive antennas, and A is the average 
Stokes 4x4 real matrix 


\m u \ 


where 

u denotes the transpose of u, 
u and are 3-element real vectors, and 
Q is a 3 X 3 real matrix, 
by maximizing the expression 

S/AS' — m + u • + v • s r + s r • Qs t 


Maximizing the signal-to-noise ratio in this case results 
in implementing the procedure developed with refer- 
ence to Equations (14) through (32) with FC 1 ) being the 
signal Stokes matrix of (e.g., a man-made target) and 
F< 2 ) being M^. A test was conducted on the man-made 
target Stokes matrix F(“) described with reference to 
Equation (33) in the previous paragraph. The optimiza- 
tion signature FIG. 8 is to be compared with the optimi- 
zation signature of FIG. 2. The effect of symmetrized 
noise is clearly noticeable in FIG. 8. The optimum po- 
larization configuration for the receive antenna and the 
transmit antenna is now i|/ r — 147.5°, 5°, 

^r= 48.9°, *r-3.8°. 

Conclusion 

The most important result of this invention is the 
development of an optimum polarization filter either for 
(1) maximization of a single target return or (2) contrast 
enhancement between two target types. Because the 
filter is based on the Stokes matrix representation, it is 
not restricted to point targets and can be applied to 
extended targets. The filter is also equally valid for 
bistatic and monostatic radar data. 

The invention is first embodied in a filter which maxi- 
mizes contrast between urban and natural targets. The 
same filter was found to be optimum when the natural 


using the Lagrangian multiplier method to produce the 
equations 

u + Qs r = Ai s t 
V + Qs { = \is r 

where X\ and X 2 are the Lagrangian multipliers, and s, 
and s r are Stokes vectors characterizing transmit and 
receive antennas, and 

having a known Stokes vector for the transmit an- 
tenna determining the receive antenna polarization 
which maximizes the signal-to-noise ratio as a third 
step from the equation 

v 4 - Qs t 

Sr ~ iiv + e*u 

where || || denotes the norm, thereby effectively maxi- 
mizing received power by matching the polarization of 
the scattered wave incident upon the receive antenna 
with polarization of the receive antenna such that 
power P is maximized under the constraints of assump- 
tions set forth in the first step above in accordance with 
the equation 


45 


50 


55 


60 


65 



13 


5 , 057,843 


14 


Pmax ~ S r T A S t - m + U • S t + V • S r + S r • Qs t . 

2. In a synthetic aperture radar imaging system as 
defined in claim 1, a method for maximizing the ratio of 5 
signal power pO) scattered by one type of target to the 
signal power p< 2 ) scattered by another type of target in 
order to maximize the ratio of the powers scattered by 
two different types of targets comprising the steps of 
first assuming the transmit antenna polarization is such 
that the Stokes vector of the transmit antenna for the 
two target types are defined as 


** ± \| ' Sfa — ~ 

N 5(a) * 5(a) 

which always represents a fully polarized wave, and 
then solving for a, an unknown, by substituting the 
expression 

±-==S=L =r* 

N 5(a) • s(a) ~ 

for s r into the equation 


ii = FV)s t « 



■Soi + s r « 51 
SQ2 + 5 r * S2 



and the Stokes vector of the receive antenna is defined 



where Sr is a completely polarized wave that satisfies 

•SfO — s r • Sr ~ 1 


to produce the quqadratic 

20 

($02 “ s 2 • S 2 )a 2 - 2(Sq\Sq 2 - s\ • 5 2 )a + (5 qj - s\ • 5i) = 0 

whereby the optimum values of a if So 2 2: t^S 2 S 2 
25 

a max — 



-501*502 - *1 * *2 

( *5 0 i5o2 ~ s\ ■ S2 

$01 — 5J • 51 

30 

+ 

$ 

1 

l *$02 ~ s 2 ’52 

J *5 02 - 52 • 5 2 


and 


and as a second step maximizing the contrast 


C(5 r ) = 


■Spl + S r - 5 ! 

5 q 2 + s r ‘ S2 


35 


a min — 


3. A method as defined in claim 2 wherein the second 
step of maximizing contrast C(s r ) subject to the con- 
straint that SfO^=SrSr— 1 is carried out by introducing 
Lagrangian multipliers and after performing differentia- 
tion finding that the optimum polarizations are the solu- 
tions to 


*Soi*So2 - Si -S2 

f *S01*S02 - 51 * 52 

| Sqi — 5i • 5i 

S 02 - s 2 * s 2 \ 

1 Sq 2 - S 2 -S2 

1 *Sq 2 - 52 • 5 2 


which define optimum values of said contrast ratio C(s r ) 
that are real, positive and definite values. 

5. A method as defined in claim 4 wherein equations 
for a max and a m i„ are substituted in equation 

Sa-as 2 = liSr 


45 


Soi + 5 r • 51 
1 S02 + Sr • 52 


= 


where /x is the Lagrangian multiplier. 

4 . A method as defined in claim 2, wherein a short- 50 
hand term 

_ j01 + s r m 5 1 

a *$02 *f S r • 52 

is introduced, and the equation 
Soi -t- s r - 5i 

51 ~~ 52 $Q2 + Sr • n ” ^ 60 

is rewritten using that shorthand term as s(a)—ixs r 
where s(a)=si — as 2 such that 

65 

fi = ± N 5(a) • 5(a) 

and the optimum polarizations are therefore 


and substituting the resulting expression for a in said 
equation 


N 5(a) • 5(a) 

to find the optimum antenna polarizations to be 

— s ^ a max) 

Smax ~ ll*a m „)ll 

and 

s( a min) 

Sm in = lKa mi »)ir 

where | | denotes the norm, and choosing the receiver 
antenna polarization to be orthogonal to the completely 
polarized scattered wave 
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*2 

s rmax — ~ v i ■ 

N Si • s 2 


for maximum contrast ratio, thereby providing an opti- 
mum receive antenna polarization for maximum con- 
trast ratio for a given transmit antenna polarization. 

10 
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6. A method as defined in claim 5 including the step 
of varying the transmit antenna polarization to maxi- 
mize contrast between targets of different types over all 
possible configurations. 

5 7. A method as defined in claim 6 wherein said differ- 

ent types of targets are man made in urban areas for one 
type and natural in nonurban areas for the other type, 
said natural type including forest, grass and ocean areas. 
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